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HEAT -TRANSFER AND TEMPERATURE DISTRIBUTION 


OF A CONICAL POROUS MEMBRANE DURING 
PLANETARY ENTRY 
By Latif M. Jiji 
Ames Research Center 


SUMMARY 


An approximate analysis is made of the stagnation point heating rate and 
surface temperature of a conical porous membrane used as a decelerator during 
planetary entry of vehicles of low ballistic parameter. Ideal gas approxima- 
tions were used. Order of magnitude variations of the porosity and ballistic 
parameter, as well as variations in cone angle, were made to find their gross 
effect on the wall temperature of the membrane. Also, the wall temperature 
was determined for an order of magnitude variation in the distance rearward 
along the membrane from the stagnation point. Porosity was found to have an 
adverse effect on the wall temperature. To avoid high wall temperatures it 
was necessary to use unrealistic values of ballistic parameter and cone angle. 


INTRODUCTION 


Decreasing the ballistic parameter, W/CpA, for a given sized vehicle 
changes the entry trajectory and decreases the maximum heating rate. A number 
of variable -drag devices have been considered for possible use during entry 
into a planetary atmosphere. In particular, deployable and variable -geometry 
configurations have received considerable attention recently. Examples are 
the^Avco Drag Brake (ref. 1), the "tension structure" (ref. 2), and the 
RotoYnet (refs. 3 an( ^ 4) . These decelerators may be competitive with the 
ablation- cooled systems when their basic structural material is lightweight, 
porous, flexible membranes. Furthermore, varying the frontal area of the 
decelerator offers a means of modulating drag for controlling the trajectory. 

The surface of decelerators is cooled by reradiation so that the 
equilibrium wall temperatures have to be within the operational level of the 
membrane material. To establish the potential value of decelerators it is 
essential to determine the aerodynamic heating and the surface temperature of 
porous membranes in a hypersonic stream. Theoretical as well as experimental 
studies of flow configurations and heat -transfer characteristics of such sys- 
tems are meager. Flow patterns and stagnation point heating of a porous mem- 
brane were investigated analytically in reference 5* A study of the heat 
transfer to a highly porous nose cap was carried out in reference 6. 



The purpose of this paper is to investigate the effect of porosity on the 
heating rate and surface temperature at the stagnation point and on the coni- 
cal skirt of a porous blunted cone during entry. Based on ideal gas approxi- 
mations, expressions are derived showing the influence of porosity and cone 
angle on the aerodynamic heating and wall temperature . Calculations are made 
to illustrate these effects for typical trajectories of vehicles with low 
W/CpA. In this report, a low value of W/CpA is taken to be 1.0 lb/ft 2 , or 
less. 


NOTATION 


A area, ft 2 

Cp drag coefficient 

C n nozzle flow coefficient 

Cp m specific heat at constant pressure, evaluated at l/2(T^. + T w ), Btu/lb m -°F 

g local gravitational acceleration, ft/sec 2 

g c gravitational constant of proportionality, 32.2 lb m -ft/lbf-sec 2 

H enthalpy, Btu/lb m 

k ratio of specific heats 

M Mach number 

m mass flow rate, lb m /sec 

p pressure, lb f /ft 2 

P porosity 

P G critical porosity 

Pr Prandtl number 

q. heat -transfer rate, Btu/sec-ft 2 

R radius of curvature, ft 

Rg gas constant, ft-lbp/lb m -°R 

R q radius of the earth, ft 

T temperature, °R 
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free-stream velocity, ft/sec 
V Q entry velocity, ft/sec 

V w suction velocity at wall, ft/sec 
VF view factor 

W vehicle weight, lb f 

X c axial distance, measured from the virtual apex of blunted cone, ft 
Y altitude, ft 

Y 0 entry altitude, ft 

3 pressure gradient parameter (3 = l/2 for a spherical stagnation point 

and 3=0 for a cone) 

7 q entry angle 

6 half-cone angle, deg 

e surface emissivity 

(i viscosity, lb^-sec/ft 2 

p mass density, lb m /ft 3 

p Q mass density at sea level, lb m /ft 3 

a Boltzmann's constant, 0.48X10 -12 Btu/sec-ft 2 °R 4 

conical skirt suction heat -transfer parameter, equation (l6) 

*i|f stagnation point suction heat-transfer parameter, equation (ll) 

s 

Subscripts 

e external to boundary layer 

f fiber 

imp impermeable 

n nose cap 

r radiation 

s stagnation point 
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sue suction (porous surface) 
t stagnation state 

w wall 

oo free stream 

* sonic section 


ANALYSIS 


The physical system under consideration is a blunted porous conical 
membrane in hypersonic flow at zero angle of attack (see fig. 1). The flow 
over a porous surface in hypersonic flight is complex; thus to analyze the 
aerodynamic heating of such a system, simplified flow patterns are assumed. 


Flow Patterns 

Flow patterns at the stagnation region of a porous nose cap vary with the 
porosity. At high porosities discrete detached shocks form over the individ- 
ual fibers as shown in figure 2(a) . As the porosity is decreased, the shocks 
interact and ultimately coalesce into a single detached shock. The critical 
porosity, P c , is defined as the minimum porosity for which there is no spill- 
age (i.e., the entire free-stream flow passes through the porous surface). 

This flow configuration is shown in figure 2(b) . At porosities lower than P c 
the flow through the pores is choked and spillage takes place as illustrated 
in figure 2(c). At some lower limit of porosity, which may be considerably 
less than P c , the detached shock must obviously approach that of a solid nose 
cap. Thus for small porosities it is possible to treat the flow near the 
surface as boundary -layer flow with surface suction. 

One-dimensional, nonviscous, compressible flow relations for an ideal gas 
are used to obtain an expression for the critical porosity, P c , at the 
stagnation point (appendix A): 


p c = 
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Equation (l) is plotted in figure 3 for various values of the ratio of 
specific heats , k , to indicate the approximate real gas effects. Since k 
decreases with increasing Mach number (up to M = 25), figure 3 indicates that 
the critical porosity for hypersonic speeds is approximately 1/3. 

The flow pattern over the conical skirt is complex; however, at low 
porosities (lower than P c ), the flow may again be treated as boundary-layer 
flow with surface suction. The detached shock over the conical surface for 
this flow pattern approaches that over the corresponding portion of a blunted 
solid cone. 


Aerodynamic Heating 

Analysis of the aerodynamic heating is based on the following two models: 

(1) Discrete shocks 

(2) Boundary-layer surface suction 

In both models the angle of attack is assumed to be zero. 

Heat transfer to fibers with discrete shocks .- In this flow pattern each 
fiber acts as a blunt body in a hypersonic stream. Since the heat -transfer 
rate is maximum at the stagnation point, attention is to be focused on this 
region of the fiber. The convective heat-transfer rate at the stagnation 
point of a sphere of radius R is given by (see refs. 7-9): 


17,000 ( Pco7' 5 

v °° 1 

SI W 

L/g(R Q + Y)-l 


3 . 15 


( 2 ) 


Using this expression limits the analysis to velocities at which the radiative 
heating can be ignored in comparison to the convective heating. 


Equation (2) may be used for comparing the heat-transfer rate at the 
stagnation point of a cylindrical fiber of radius Rf with that of a solid 
nose cap of radius R n . For a cylinder and a sphere of the same radius, 
equation (2) differs only by a factor of l//2". The resulting expression for 
the ratio of the fiber to nose cap heat transfer is: 



<ln 



(3) 


For the range of nose cap to fiber radii considered, the heat-transfer ratio 
can become very large, as shown in figure 4. It is clear that the heat- 
transfer rate to the individual fibers is much greater than that at the stag- 
nation point of a nose cap. Therefore, it becomes obvious that porosities 
larger than the critical value cannot be tolerated. 
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Heat transfer to porous surface with boundary-layer suction. - For 
porosities less than the critical porosity, P c , the flow oyer a porous surface 
may be approximated by bound ary-layer flow with surface suction. Limited work 
on this problem has been done for supersonic flow. A theoretical analysis of 
the effect of mass injection through the boundary layer on heat .transfer 
around blunt bodies of revolution is given in reference 10. In reference 11 
the analysis was applied to suction by changing the sign of the mass flow 
term. This theory was found to agree well with experimental findings 
(ref. 11). For laminar flow of partially dissociated air, the effect of suc- 
tion on heat transfer, as verified in reference 11, is given by: 




L suc _ 


= 1 + 


^lnrp 


vhere 


or 


p w v. 

q. 


B = OLJL (H t - H w ) 


imp 


B ' A 2 c pJ T t - T »> 

imp 


00 


(5) 


Equation (4) was correlated in reference 11 for values of the suction 

par ame ter B up to 4.4. Since for porous surfaces it is not known a priori 

that T w /T t « 1, the wall temperature has been retained in equation (5). The 
parameter (3 takes into consideration the effect of pressure gradient outside 

the boundary layer. The stagnation temperature T t is given by 


T 


t 




(6) 


The suction mass flow rate per unit area, p w V w , is obtained by approximating 
the flow through the pores with one -dimensional, compressible nozzle flow and 
assuming that the pressure downstream is low enough for conditions to be sonic 
at the throat. When the continuity equation is applied at the wall (i.e., the 
inlet section of the nozzle) 

= PwAAw 


or 


_SL = n V ^2- 
a Pw w A 
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Noting that the porosity P is defined as 


P 





the above becomes 


o v » JL p 
M w w 1 * ' 


( 7 ) 


The maximum flow per unit area, m/A^, is given by 


k+i 

m - c &£ fl_2_ f* 

a ; - c *Ir s [k + i ) ^ 


( 8 ) 


Equations (4), (5)> (7)> and (8) give 


/ k+x 

- - 1 + W" c i^^; )Pt " (^i^) 


i inrp 


(9) 


In equation (9) Pt w is * t5ie stagnation pressure on the porous surface. This 

may be approximated by the stagnation pressure downstream of the detached 
shock. 


Equation (9) is now applied to two regions: (a) the stagnation point and 

(b) the conical skirt. Expressions for p^ and q^^p must be obtained. 


Stagnation point: The detached shock at the stagnation region of the 

blunted nose cap may be approximated by a normal shock. The stagnation pres- 
sure downstream of the shock, p+- , is given by 

w k 



( 10 ) 


The heat transfer to an impermeable surface at the stagnation point is 
given by equation (2). 

To determine the effect of suction on the stagnation-point heat transfer, 
the quantity is defined as 
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( 11 ) 
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Equations (2), (9), and (10) give 


i s = 0.362X10' 4 yp. 
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( 12 ) 


In equation (12) the Prandtl number for air is taken as equal to 0.71 and 0 
at the stagnation point is 0*50. The dimensionless ratio (p e [i e /p^ v )°'° 4 is 
assumed to be unity. This is justified since (PeM-e/Pw^w) is orcier of 
magnitude of one. 


Equation (12) gives the value of as a function of altitude for a 

vehicle of known trajectory. 


Conical skirt: Since streamlines are approximately parallel to the 

conical surface, equation (8) overestimates the suction mass flow rate if the 
actual stagnation pressure, pt w > is used. A more realistic estimate can be 
obtained by replacing Pt w in equations (8) and ( 9 ) by the static pressure at 
the wall, p v . If the effect of blunting is neglected, then from Newtonian 
flow approximations the static pressure over the conical skirt is given by 


p 

— = 1 + kJVr sin 2 6 

p 00 

00 


( 13 ) 


For cones in dissociated air the shock may remain attached for half angles up 
to 6o°. The Newtonian approximation underestimates the actual value of p w 
by approximately 6 percent for values of 5 ranging from 0° to 6o° (ref. 12). 

The heat-transfer rate to an impermeable cone is correlated with the 
free-stream velocity, density, and cone angle in reference 13 . The following 
correlation equation applies to laminar flow over pointed cones at zero angle 
of attack: 
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Since blunting has a minor effect on the heat-transfer characteristics when 
the distance X c is measured from the virtual apex of the blunted cone 
(ref. 7)* equation (l4) can be applied to blunted cones. 

When equations (13) and (l4) are substituted into equation ( 9 ), and since 
0 = 0 for a cone, the following is obtained: 


0 5 / k+1 

= 1 . 08 x 10 “ 3 ( Bsa'l ( 421^ / g c k ( 2 j 1 C P m ^ Tt “ T 0 (f + sin 2 6)p m 
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sue 


C n P n/x^ ^ q imp 


(15) 


(16) 


Equilibrium Wall Temperature 

To determine the equilibrium wall temperature, the convective heat rate 
is equated to the heat loss by radiation* Radiation heat- transfer rate is 
given by 


q r = e(VF)arJ (17) 

When equation (2) is written for flow over a cylinder and is combined with 
equation ( 17 ), for discrete shocks. 


L w 
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17,000 ( Poo 
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( 18 ) 


For flow with boundary layer suction the wall tenqoerature at the stagnation 
point of a porous nose cap is given by 


17,000 
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Equation (19) is obtained from equations (2), (11), (12), and (17)* Similarly, 
for the conical skirt, equations (l4), (15), (l6), and (17) give: 


222 


/ \ °* 5 / N 3 

(0 &) 


sin 25 + 


0.576 C n P 



f" 1 Cp m (Tt - T w ) Poc 


■JT t 


Q 2 

(1 + sin 5 ) = e(VF)a'I* 


( 20 ) 


APPLICATIONS 


To investigate the effect of porosity on the heat-transfer rate and sur- 
face temperature distribution, the analysis is applied to vehicles with low 
ballistic parameters. Three trajectories based on the 1959 ABDC atmosphere 
(ref. 14) were computed for values of w/CpjA of 0.01, 0.1, and 1.0 lbf/ft 2 . 
Figure 5 shows these trajectories for an equatorial entry at 25,000 ft/sec, 
500,000-ft initial altitude, and with an initial entry angle of -5°. 

The suction heat-transfer parameter at the stagnation point of a 

porous nose cap determined from equation (12) is plotted in figure 6 for a 
wall temperature of 1000° R. For the conical skirt, equation (16) is plotted 
in figure 7 for half -cone angles of 15°, 45°, and 6o°. Note that the ballis- 
tic parameter W/OqA has an appreciable effect on both tyg and ty c . 

Stagnation wall temperature profiles for fibers with individual shocks 
are plotted in figure 8 (eq. (l8) ) for fiber radii of 0.0025, 0.0125, and 
0.025 inch. Even for a ballistic parameter as low as 0.01 lbf/ft 2 the result- 
ing temperatures are excessively high. Thus configurations which result in 
discrete shocks over the individual fibers cannot be tolerated. 

For a porous nose cap with a boundary-layer-suction flow pattern, 
temperature profiles at the stagnation region of a 2-foot radius nose cap are 
plotted in figure 9* Profiles are shown for values of C n P ranging from 0 
to 0.60; C n P = 0 corresponds to a solid surface, while C n P = 0. 60 is an 
approximate upper limit corresponding to the critical porosity for k = 1.4. 
For values of C n P = 0.01, the effect of suction on the temperature distribu- 
tion is negligible. However, for C n P =0.1 the effect is significant, par- 
ticularly for a ballistic parameter of 1.0 lbf/ft 2 or greater. It is 
therefore desirable to have as inefficient a nozzle as possible, even for low 
values of porosity. 

Figure 10 gives the effect of ballistic parameter on the maximum 
stagnation-point wall temperature for various values of C n P* The wall tem- 
perature is observed to increase rapidly with w/lpA. For reasonable values 
of W/Cj)A for this concept, about 1.0 lbf/ft 2 , C n P of 0.1 increases the max- 
imum wall temperature relative to the nonporous value over 50 percent. On the 
other hand, for w/C-qA = 0.01 lbf/ft 2 the corresponding increase in the 
maximum temperature is about 10 percent. 
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Temperature profiles at the conical skirt for X c = 10 feet are shown in 
figure 11 for various values of C n P, w/Cj)A, and 8. At 8 = 6o° (fig. 11(c)) 
the porosity effect is significant for c n p as low as 0.10 for 
W/CpA = 1.0 lb f /ft 2 . These profiles indicate that the adverse effect of 
porosity on the temperature becomes more severe as w/CpA and/or 8 are 
increased. 

Figure 12 shows the variation of maximum wall temperature with X c . It 
is observed that the effect of changing X c on the wall temperature is less 
pronounced when the effect of suction is appreciable (i.e., at high values of 
W/CpA and/or 8). 

Temperature distribution charts for both the nose cap and the conical 
skirt indicate that if W/GqA could be reduced to 0.01 lbf/ft 2 , the maximum 
wall temperature could be maintained at reasonable low levels. However, in 
this case the ratio of payload to weight for this type of retardation loses 
its advantage over that of the ablative heat shield with parachute recovery 
system. 

A major factor in the design of a rotating conical skirt is the size 
required. For example, to recover a payload of 500 lb (including decelerator 
weight) with C;n = 1 would require a skirt 25 feet in diameter for 

W/C pA = 1.0 lbf/ft 2 . If, to avoid the severe heating problem, W/CpA were 

reduced to 0.01 lbf/ft 2 , the diameter would increase to 250 feet. These cal- 
culations were made for a value of k of 1.4. While real gas effects can be 

evaluated only from real gas properties, rather than ideal gas relations, lim- 

ited calculations were made for k = 1.1 in an attempt to estimate real gas 
effects. Wall temperatures were found to be approximately 15 percent lower 
for k = 1.1 than for 1.4. 


CONCLUDING REMARKS 


The wall temperature for a porous conical membrane during planetary entry 
is excessively high when discrete shocks form over the individual fibers; thus 
membrane porosities must be less than the critical porosity. For the bound- 
ary layer suction flow pattern, associated with low porosities, and reasonable 
values of W/CpA (about 1.0 lbf/ft 2 ) the temperature increase due to porosity 
is too severe even for relatively low porosity. For lower values of w/C-qA, 
the payload-to-weight ratio loses its advantage over that of the ablative heat 
shield with parachute recovery. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, Feb. 7, 1967 
124-07-02-27-00-21 
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APPENDIX A 


CRITICAL POROSITY 


For one -dimensional isentropic 
unit area is given by: 


flow of a perfect gas, the mass flow per 



A 



?t 

-Th 




k+i 

\a(k-i) 


(Al) 


If equation (Al) is applied to the flow upstream of the shock, 


m = 

A w 


^c P t. 





1 + M 2 

2 00 


k+i 

a(k-i) 
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If it is assumed that the pressure downstream of the porous surface is low 
enough to give sonic speed at the throat, equation (Al) gives 
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m _ / *S C ( 2 Y' 1 *t 
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Equations (A2) and (A3) give the critical porosity 
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(Alt) 


The stagnation temperature ratio T # /T t is equal to unity. The stagnation 

"t 00 

pressure ratio is obtained from normal shock relations 
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Equations (Ah) and (A5) give 


l 
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Figure 1.- Conical membrane in hypersonic flight. 
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(a) 5 = 15° 

Figure 7*“ Heat -transfer parameter for conical skirt with suction, 


Yo 

Vo 

Yo 

T w 


= 500,000 ft 
= 25,000 ft/se. 
= -5° 

= 1000 °R 



(b) 5 = ^5° 


Figure 7 .- Continued. 












Maximum wall temperature, 



w 

CdA' 


Ibf/ft 2 


Figure 10.- Maximum stagnation point wall temperature; e (VF) = 1.5, 

R n = 2 feet. 
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Figure 12.- Effect of X c on the maximum wall temperature; € (VF) = 1.5 

C n P = 0 . 1 . 
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